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ABSTRACT 

New optical spectra of 28 Hll regions in the MlOl disk have been obtained, yielding 10 new detections 
of the [Olll]A4363 auroral line. The oxygen abundance gradient measured from these data, combined with 
previous observations, displays a local scatter of 0.15 ±0.03 dex along an arc in the west side of the galaxy, 
compared with a smaller scatter of 0.08 ±0.01 dex in the rest of the disk. One of the H II regions in our sample 
(H27) has a significantly lower oxygen abundance than surrounding nebulae at a similar galactocentric distance, 
while an additional, relatively nearby one (HI 28) was already known to have a high oxygen abundance for its 
position in the galaxy. These results represent marginal evidence for the existence of moderate deviations from 
chemical abundance homogeneity in the interstellar medium of MlOl . Using a variety of strong-line abundance 
indicators, we find no evidence for significant large-scale azimuthal variations of the oxygen abundance across 
the whole disk of the galaxy. 

Subject headings: galaxies: abundances — galaxies: ISM — ISM: abundances — galaxies: individual (MlOl) 



1. INTRODUCTION 

Radial gradients of metallicity in spiral galaxies have been 
known for a long time, with typical values of 0.03-0.10 dex 
kpc~' . The study of giant extragalactic H II regions offered 
the first evidence for the prese nce of oxygen abundance gra- 
dients in galaxies ( Searle 1971| l and has been extensively used 
to probe the stellar populations an d the chemical composi- 
tion of nearby star-fo rming galaxies (Vila-Costas & Edmunds 
T992l |Zaritsk y et al.'1994; van Zee et al. 1998; Bre solin et al. 



1999 I, placing strong constraints on galactic chemical evolu- 



tion models (Chiappini et al.^2001 ,,Fu et al.|2009) . In contrast, 
surprisingly Uttle is known about the possible presence of az- 
imuthal asymmetries in the abundance distribution in spiral 
disks. From a theoretical point of view, a virtually uniform 
distribution, with negligible scatter, could be expected along 
the azimuthal direction, as a resul t of relatively fast mixing 
processes (except near corotation, [Scarano & Lepine]|2013| l, 
with timescales on t he order of ^1 00 Myr, in the turbulent 
interstellar medium (IRoy & Kunth 1995; Yang & Krumholz 
[20l2| . 

Deviations of the radial abundance gradient from a sim- 
ple exponential seem to be well -established from stu dies of 
Cepheids (e.g. |Luck et al.| 2003| l and open clusters ( Twarog| 
et al.|199T Lepine et al.|201 1 1, showing a discontinuity at a 



galactocentric distance of ~8.5 kpc and a flattening beyond 
that radius, possibly caused by the barrier effect of corota- 
tion, which isolates the inner and outer regions of the disk 
one from the other , due to opposite directions of gas flow 
( Lepine et al.|2011] l. Similar features can be observed in the 
H II region abundance distribution of nearby spira l galaxies 
( [Bresolin et al. 2009b 2012; Sca i'ano et al.|[20n] l. Investi- 
gations of Cepheid s (Pedicelli et al. 2009| l and Hll regions 
( Balser et al.pOll) in the Milky Way have provided indica- 
tions for the presence, at least near the corotation radius, of 
azimuthal gradients (^0.05 dex kpc"') and chemical inhomo- 
geneities, which can be attributed to the spiral arm structure 
and the consequent non-uniformity of the spatial distribution 
of gas and star formation in the disk of the Galaxy (Lepine] 
et al. 2011| . These results stress the importance of includ- 



ing the azimuthal coordinate, in addition to the radial one, in 
modeling the chemical evolution of the Milky Way. 

The situation in other spiral galaxies is less clear, and suf- 
fers from poor statistics (small number of H II regions stud- 
ied in a give n galaxy) and mea surement uncertainties. In the 
galaxy M33, |Rosolowsky & Simo n (2008) measured a sub- 
stantial intrinsic dispersion of OJTT dex (in addition to the 
dispersion due to observing errors) in the Hll region oxy- 
gen abundance at constant radius. However, Bresolin (2011) 
showed that this is due to poor detections of the [O lll]A4363 
auroral line used to derive nebular electron temperatures, and 
that a much smaller scatter of ^0.06 dex, consistent with the 
measurement uncertainties, is obtained from the best-quality 
data. In a wide-field integral field spectroscopy study of 



NGC 628, |Rosales-Ortega et aLl ( [20Tl] found that the radial 
metallicity gradient varies slightly for different quadrants, al- 
though the differences are comparable to the uncertainties in- 
troduced by the nebular abundance diagnostics and can be re- 
lated instead to systematic variations of the ionization param- 
eter. The narrow-band imaging across the disk of the same 
galaxy by Cedres et al. (2012 1 suggests that the detection 
of chemical inhomogeneities may depend on the choice of 
strong -line diagnostics used to derive the oxygen abundances. 
Recently, a large degree of inhomogeneity on relatively small 
spatial scales (^0.5 kpc) has been claimed by Sanders etaD 
(20121 for the Hll region oxygen abundances in M31. ITiis 
result is not confirmed by other investigations in the same 



galaxy, including the recent one by Zurita & BresoUn ( [2012 1, 
and it is in general still unclear what role the sample selection 
and choice of chemical abundance diagnostics play in these 
results. 

In an attempt to provide new observational constraints 
on the presence of chemical inhomogeneities and azimuthal 
gradients in nearby spiral galaxies, we turned to MlO l 
(NGC 5457), a nearby (D = 6.85 Mpc, |F;reedman et al.|2001|), 
nearly face-on, grand-design spiral galaxy, that has been 
widely studied as a prototype system for the investigation of 
radial gra dients of eleme nt abundances (see Kennicutt et al. 
[2003 1 and [Bresolin|2007| for a review and references). ITiere 
are more than 10^ cataloged H ii regions in its disk ( Hodge 
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et al.|1990| l. [Kennicutt & Garnett| ( |1996| suggested the pres- 
ence of a possible azimuthal asymmetry in the oxygen abun- 
dance distribution between the southeast (SE) and northwest 
(NW) regions of MlOl, based on strong-line abundance di- 
agnostics. However, this result is subject to the uncertainties 
and dis crepancies concerning the calibrations of t hese diag- 
nostics dBresoUn||2007t [Kewley & Enison|[2008l ). I Kennicut t 
et al. (2003D improved on previous studies usmg direct oxy 



gen measurements for 20 H II regions, which rely on the deter- 
mination of the electron temperature from auroral-to-nebular 
line ratios, such as [O lll]A4363/A5007. These authors found 
a ^0.2 dex metallicity spread for H II regions located at a sim- 
ilar galactocentric distance in the southwestern section of the 
galaxy. However, limited by the size of their H II region sam- 
ple, they could not confirm the large-scale azimuthal asym- 
metry between the SE and NW, suggested by Kennicutt &| 
|Garnett'(T996^. 

In this paper, we present new spectroscopic observations of 
28 H II regions in MlOl, and are able to derive direct oxygen 
abundances for a subset of 10. This new dataset allows us to 
address the possibility of a non-axisymmetric distribution of 
the nebular oxygen abundance in this galaxy. In addition, be- 
cause MlOl is known to be currently experiencing an infall of 
high-velocity gas, as shown by H I maps (Sancisi et al. 2008 ), 
and has likely been recently subjected to interaction events, 
as shown by the lopsidedness of i ts disk, the peculiar spiral 
structure and faint tidal structures (! Waller et al.|1997"{ |Mihos| 
|et al. 2012), it is well suited to verify to what extent these 
events can affect the distribution of metals in a spiral disk. 

This paper is organized as follows. In Sectionl2]we describe 
the observations and data reduction procedureslln Section |3] 
we present and analyze the radial gradients, and the local az- 
imuthal variations of oxygen abundance, using both the direct 
method ( [O III] A4363-based) and strong-line diagnostics. We 
test the asymmetry of the oxygen abundance distribution be- 
tween the east and west sections of MlOl in Section |4] and 
present our conclusions in Section|5] 

2. OBSERVATIONS 

We obtained multi-object spectroscopy of MlOl during the 
nights of May 7-8, 2010 with the R-C spectrograph at the 
Mayall 4-m telescope of the Kitt Peak National Observatory, 
using multi-object masks in five different fields, with 2 arcsec- 
wide slits. The data were acquired at airmasses smaller than 
1.25, to minimize the impact of differential atmospheric dis- 
persion. In each field we integrated for 3 x 2400 s using the 
KPC-lOA grating (2.75 A pixel"'), which yielded a --7.5 A 
FWHM spectral resolution and covering approximately the 
wavelength range between 3600 A and 7500 A. 

Standard IRAlQ tasks, in combination with the PyRApF] 
command language, were used for bias subtraction, flat-field 
correction, cosmic ray removal, spectral extraction, image co- 
addition and wavelength calibration (using He-nNe-nAr lamp 
frames). We extracted one-dimensional spectra for each tar- 
get Hll region after executing all the standard reduction 
procedures. Three standard stars (Feige 34, Feige 67 and 
BDh-284211) were observed several times each night for the 
flux calibration. 

' IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 

^ PyRAF is a product of the Space Telescope Science Institute, which is 
operated by AURA for NASA 



Our observations yielded spectra for 28 H II regions. Their 
celestial coordinates, galactocentric distances normalized to 
the isophotal radius Rq, position angles (f> relative to the galac- 
tic center and identifications from the Hodge et al. ( 1990) cat- 
alog, are presented in Table [1] (where objects are listed in or- 
der of increasing declination). We have adopted the follow- 
ing parameters: center coordinates RA = 14:03:12.5, Dec = 
-1-54:20:56 (J2000), position angle of the major axis = 37 deg, 
inclination angle =18 de g (de Blok et al. 2008 ), distance = 
6.85 Mpc (Freedman et al. 2001 ) and isophotal radius TJq = 
14.'4 (de Vaucouleurs et al. 1991), corresponding to 28.7 kpc. 
Fig. [T] shows the location of the target Hll regions in MlOl. 
Hll regions for which we detected the [Oni]A4363 auroral 
line are labeled by double symbols, and we identify the differ- 
ent data sources considered in this paper with different sym- 
bols. 

2.1. Line flux measurement 

The emission Une intensities were measured with the splot 
routine in IRAF, by integrating the flux under the line pro- 
files between two continuum points selected interactively. A 
multi-gaussian profile fit was performed if evident line blend- 
ing occurred. 

A correcti on for interstell ar reddening was performed 
adopting the Howarth (1983i analytical formulation of the 



Seaton ( 1979[ ) law, assuming a total-to-selective extinction ra- 
tio Ry = AvIEb-v = 3.1. By comparing the measured inten- 
sities of Ha and H7, relative to H/3, to the c ase B theoretical 
values taken from |Storey & Hummer| ( [T995| l calculated at the 
electron temperatures determined from the auroral lines (or 
10,000 K if these lines were not detected), we obtained the 
reddening coefficient c(H/3). The Balmer lines ratios (Ha/H/3, 
H7/H/3) were corrected for underlying stellar absorption in an 
iterative manner. 

The resulting reddening-corrected line flux measurements, 
normalized to H/3 = 100, are tabulated in Table |2] The errors 
in these fluxes were estimated from the uncertainties in the 
line intensity measurements and the flux calibration, the scat- 
ter of the continuum near the emission lines, and the uncer- 
tainties in the extinction coefficient. We note that, in general, 
the H II regions in our sample are extended, and are not fully 
covered by the width of our slits. Thus, the H/3 fluxes reported 
in Table 2 should be considered lower limits to the real fluxes. 

2.2. Electron temperatures and oxygen abundances 

Electron temperatures and oxygen abundances were ob- 
tained adopting a two excitation zone structure for the H II re- 
gions, characterized by the temperatures T [O ll] and T [O III]. 

To calculate the total oxygen abundances, we made the 
usual assumption: 



0/H = (0-' + 0-'+)/H+ 



(1) 



The electron temperature {T^) for the 0''"^-emitting re- 
gion (T[Olll]) and for the 0"^-emitting region (T[Oll]) 
were derived from the auroral-to-nebular line intensity ratios 
( [Oiii]A4363/AA4959, 5007 and [ Oii]A7325/A3727), wh en 
available, with the nebular module ( jShaw & Dufour|1995| l of 
the STSDAS package running under IRAF. For nebulae with- 
out auroral line detections we adopted T^, = lO'' K. One or both 
of the [Olll]A4363 and [Oll]A7325 auroral lines could be 
measured in 10 out of the 28 H II regions in our sample: seven 
objects with both lines detected, three with only [O lll]A4363 
detected and one with only [Oii]A7325 detected. We com- 
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Fig. 1 . — Location of the H II regions in our study on a narrow-band Ho image of MlOl (from the Canada-France-Hawaii Telescope data archive). The whit e 
cross marks the center of the galaxy. The symbols indicate th e different data sources: blue triangl es - new ob servations; green circles -IKennicu tt et al.|j2003) ; 
magenta stars - Bresolin 1 2007 1; orange diamonds - [Kennicutt & Garnett ( 1996 1; white squares - v an Zee et al.j l 1998 1. H II regions with [O IIIJA4363 detections 
are labeled by doubled line symbols. Note that four distant objects (H1262, H681, H682 and SDHi'li) are outside the field of view (the arrows indicate the 
directions to these nebulae). The two curved dashed lines in the west mark the position of 'arc A'. 



pared our measured T [O ll] values with those obtained from 
the theoretical relation obtained by Garnett ( 1992): 



T[0 II] = 0.70 T[0 III] -I- 3000 K 



(2) 



and found that the two were consistent with each other, within 
the uncertainties. However, since the T [O n] values ob- 
tained from the auroral-to-nebular line ratios have signifi- 
cantly larger errors (between 300 K and 900 K) than those 
propagated from T [O III] through Eq. |2l( 170 K to 400 K), a 
result already noticed by Kennicutt et aLl(!2003||, we decided 
to adopt the T [O ll] values calculated via Eq. [2] In the end, 
this choice has a limited impact on the abundance determi- 
nation, since the O/H ratios obtained by adopting one or the 
other set of T [O n] values differ only by a few hundreds of a 
dex (0.07 dex at most), except for H59, where the difference 



is 0.13 dex. 

The measured values of the O'^-emitting region are 
presented in Table |2] We assume, as indicated by the 
[S II] A6717/A6731 line ratios, that aU H II regions are in the 
low density limit (n^ < 10^ cm"^). 

We derived the oxygen abundances both from the mea- 
surement of the electron temperature ('direct' method or Tg 
method), obtained as explained above, and strong-line diag- 
nostics. Strong-line diagnostics were put forward to obtain 
the abundances in H II regions when the temperature-sensitive 
lines are not accessible (e.g. in metal-rich regions or high red- 
shift galaxies), by using combinations of more easily mea- 
sured strong emission lines as indicators of the abundance, 
such as the popular [O ll] A3727+ [O lll]AA4959, 5() 07)/H/3 
line ratio, otherwise known as R23 ( jPagel et al.|1979| l. As is 
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TABLE 1 
Observed H II region sample 



ID-' 


R.A.'' 


Decl." 




0' 


Other id' 




(J2000.0) 


(J2000.0) 


(degree) 




1 


14 02 46.92 


54 14 50.1 


0.50 


212 


H219 


2 


14 02 53.68 


54 15 22.8 


0.43 


206 


H321 


3 


14 02 57.54 


54 15 59.8 


0.37 


204 


H370 


4 


14 02 30.61 


54 16 09.7 


0.54 


232 


H149 


5 


14 02 25.73 


54 17 34.2 


0.53 


244 


H103 


6 


14 02 49.83 


54 17 43.0 


0.32 


226 


H260 


7 


14 02 48.78 


54 17 56.3 


0.32 


229 


H237 


8 


14 02 18.22 


54 18 36.6 


0.58 


254 


H41 


9 


14 02 27.01 


54 19 33.0 


0.48 


258 


H120 


10 


14 02 29.49 


54 19 47.7 


0.45 


260 


HMO 


11 


14 02 19.92 


54 19 57.1 


0.55 


263 


H59 


12 


14 02 16.85 


54 21 00.1 


0.58 


271 


H27 


13 


14 02 22.06 


54 21 57.7 


0.54 


278 


H79 


14 


14 02 17.83 


54 22 32.4 


0.59 


282 


H378 


15 


14 02 17.83 


54 22 32.4 


0.59 


282 


H37S 


16 


14 02 20.39 


54 23 13.7 


0.58 


287 


H68 


17 


1 4 02 27 50 


54 27 OS n 


0.66 


313 


H125 


18 


14 03 50.86 


54 27 35.9 


0.61 


40 


H1151 


19 


14 03 49.98 


54 27 50.5 


0.61 


38 


HI 148= 


20 


14 03 49.98 


54 27 50.5 


0.61 


38 


H1148S 


21 


14 03 49.18 


54 28 10.0 


0.63 


36 


HI 146 


22 


14 04 49.80 


54 28 14.0 


1.12 


63 


H1262 


23 


14 03 46.84 


54 28 37.3 


0.64 


33 


H1137 


24 


14 02 39.30 


54 29 04.9 


0.69 


329 


H188 


25 


14 02 49.65 


54 29 07.0 


0.64 


338 


H253 


26 


14 03 44.30 


54 29 14.4 


0.66 


29 


H1121 


27 


14 03 43.26 


54 29 49.4 


0.69 


27 


H1118 


28 


14 03 13.70 


54 35 47.0 


1.05 


1 


H681 



H II region identification (in order of increasing declination) 
^ Rigiit ascension. Units are iiours, minutes and seconds 
Declination. Units are degrees, arcminutes and arcseconds 

Deprojected galactrocentric distance in units of the isophotal radius, Rq = 14' 4, taken 
from de Vaucouleurs et al. i 1991 i 



' Position angle relative to the galaxy ce nter, in degrees 
Identification from Hodge e t al.|jl99oj 
' Two objects in the same slit 



well known, strong-line indicators, like R23, can suffer from 
degeneracy when calculating the chemical abundances, and 
abundances derived from different strong-line methods and/or 
calibr ations can vary system atically up to Alog(0/H) = 0.7 
dex ( jKewley & Enison|2008| l. 

2.3. The enlarged sample 

In order to increase the sample size and to cover a larger 
radial and azimuthal range, we also collected data from addi- 
tional sources: 



• [vanZee etal.| ( | 19981 = vZ98) published spectroscopic data 
for 13 H II regions in MlOl. These authors did not follow the 



• Kennicutt & Garnett ( 1996 =KG96) presented spectro- 
scopic data for 41 H II regions, including some distinct knots 
in large complexes (NGC 5462, NGC 5447 and NGC 5471), 
which could be used as control samples to examine the intrin- 
sic metallicity variation, if we assume that within their volume 
the oxygen abundance is homogeneous. 



object identification from the catalog of |Hodge et al. ( 19901. 
Based on their map of slit positions, w e identified the corre- 
sponding objects in |Hodge et al.| ( [T990| . 

In assembling our catalog of line fluxes for H II regions in 
MlOl, we adopted the more accurate line fluxes given by 
KBG03 and B07 for the 26 objects in common with KG96. 
Similarly, for H67 and H188, which were observed by vZ98, 
we adopted our new data or those by KBG03. 

Since in KG96 and vZ98 only the sum of the 
[Olll]AA4959,5007 and [N ll]AA6548,6583 line fluxes are 
given, we calculated the individual line fluxes of the Unes in 



• [Kennicutt et~aL] ( |2003l = KBG03) observed 26 H II re- 
gions, with 23 targets in common with KG96; 17 of those 
23 observations yielded [O lll]A4363 measurements. In ad- 
dition, they also observed 3 new objects: H70, H71 and 
SDH323 (following the designation in KBG03), all of them 
with [O lll]A4363 measurements. 



• |Bresolin| ( p007] =B07) observed 4 Hll regions (H1013, 
H493, H507 anarH972) in the inner (central 3'), metal-rich 
zone of MlOl; two (H1013 and H493) provided reliable 7; 
measurements. 



the doublets using the theoretical ratios (Storey & Zeippen 
2000). 

Table [3]presents line fluxes from different authors (KBG03 
and B07) for the six objects in common with our new dataset 
of 28 H 11 regions (labeled as "New" in the table). In Fig. |2] 
we display a comparison of line fluxes and derived quantities 
(R23, O/H). 

The agreement between the different studies is generally 
satisfactory. There is only one object, H140, with a consid- 
erable difference (~20%) in the [S ll] flux between our new 
observations and those by KBG03. Fig.|2|c) indicates that R23 
values are consistent between the different studies, with dif- 
ferences < 0.12 dex. There are only three objects in common 
that also have an [O lll]A4363-based oxygen abundances. For 
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TABLE 2 

Dereddened line fluxes, temperatures and oxygen abundances from the direct method 



ID 


Name 


[O II] 


[O III] 


[0 111] 


[Nil] 


[S II] 


F(H/3) 




c(H/3) 


r/~\ tit! 
1 [U IIIJ 


I2+IOg(W/ri; 






A3727 


A4363 


A5007 


AO JO J 


\ \ *i'7 1 T ATI 1 

AAo/ 1 /,0/iI 


(erg s^^ cm 


) 


(mag) 


(K) 




1 


H219 


413 ± 21 


0.91 ± 0.05 


178 ± 8 


36.0 ± 1.5 


41.3 ± 1.3 


3.8 xlO- 


14 


0.31 


9420 ± 258 


8.42 ± 0.08 


2 


H321 


264 ± 13 




237 ± 11 


29.5 ± 1.2 


24.7 ± 0.7 


3.3 X ur 


14 


0.43 






3 


H370 


283 ± 14 




129 ±6 


78.3 ± 2.8 


54.4 ± 1.4 


1.6 X 10- 


14 


0.53 






4 


H149 


238 ± 12 


1.68 ± 0.08 


329 ± 15 


35.6 ± 1.5 


39.4 ± 1.2 


1.6 X 10" 


13 


0.35 


9411 ± 247 


8.40 ± 0.08 


5 


H103 


343 ± 17 




158 ± 7 


42.6 ± 2.0 


38.6 ± 1.3 


1.3 xlO- 


14 


0.18 






6 


H260 


245 ± 12 




207 ±9 


36.6 ± 1.4 


35.4 ± 1.0 


3.6 X 10- 


14 


0.45 






7 


H237 


193 ± 10 




106 ±5 


47.5 ± 1.7 


38.5 ± 1.0 


8.2 X 10- 


14 


0.51 






8 


H41 


169 ±9 


3.92 ±0.18 


476 ± 21 


8.4 ± 0.4 


15.8 ± 0.5 


3.7 xlO- 


14 


0.28 


10894 ± 327 


8.25 ± 0.08 


9 


H120 


419 ±21 




82 ±4 


63.3 ± 2.8 


86.9 ± 2.8 


1.9 xlO- 


14 


0.26 






10 


HMO 


352 ± 18 




145 ±7 


57.0 ± 2.5 


67.1 ± 2.1 


3.1 xlO- 


14 


0.29 






11 


H59 


305 ± 15 


2.30 ±0.11 


290 ± 13 


17.2 ±0.8 


25.9 ± 0.9 


2.1 xlO- 


14 


0.16 


10769 ± 322 


8.24 ± 0.08 


12 


H27 


341 ± 17 


3.64 ±0.17 


251 ± 11 


21.6 ± 1.0 


58.0 ±2.0 


2.1 xlO- 


14 


0.16 


13401 ± 499 


7.96 ± 0.08 


13 


H79 


226 ± 1 1 


3.04 ±0.15 


330 ± 15 


15.5 ± 0.7 


27.1 ± 0.9 


1.7 xlO- 


14 


0.25 


11313 ± 563 


8.14± 0.11 


14 


H37 


236 ± 12 




333 ± 15 


19.2 ±0.9 


22.1 ± 0.8 


8.6 X 10- 


15 


0.19 






15 


H37 


234 ± 12 


3.31 ± 0.16 


387 ± 17 


14.9 ± 0.6 


28.5 ± 0.8 


2.5 X 10- 


14 


0.38 


11043 ± 344 


8.22 ± 0.08 


16 


H68 


391 ± 19 




171 ± 7 


33.8 ± 1.5 


75.6 ± 2.5 


1.1 xlO- 


14 


0.26 






17 


H125 


350 ± 18 




327 ± 15 


29.0 ± 1.4 


33.3 ± 1.3 


1.3 xlO- 


15 


0.15 






18 


H1151 


397 ± 21 




158 ±7 


25.0 ± 1.0 


35.4 ± 1.1 


2.1 xlO- 


15 


0.43 






19 


H1148 


177 ± 9 




293 ± 13 


17.5 ± 0.8 


25.6 ± 0.8 


7.9 X 10- 


15 


0.28 






20 


H1148 


199 ± 10 




383 ± 17 


11.0 ± 0.4 


14.5 ± 0.4 


1.5 xK)- 


14 


0.40 






21 


H1146 


212± 11 


3.07 ± 0.15 


457 ± 20 


17.7 ±0.8 


27.1 ± 0.9 


9.6 X 10- 


15 


0.27 


10219 ± 299 


8.35 ± 0.08 


22 


HI 262 


200 ± 10 




316± 14 


13.0 ± 0.7 


31.6± 1.1 


3.4 X 10- 


15 


0.22 






23 


H1137 


378 ± 19 




25 ± 1 


43.7 ± 2.0 


87.5 ± 2.9 


2.3 xlO- 


15 


0.25 






24 


H188 


682 ± 34 




164 ±7 


23.8 ± 1.1 


58.7 ± 2.0 


3.9 X 10- 


15 


0.23 






25 


H253 


430 ± 22 




73 ± 3 


40.8 ± 1.9 


84.6 ± 2.9 


3.6 X 10- 


15 


0.17 






26 


H1121 


336 ± 17 




217 ± 10 


29.3 ± 1.3 


47.7 ± 1.6 


5.7 xlO- 


15 


0.26 






27 


H1118 


222 ± 1 1 


5.97 ± 0.29 


559 ± 25 


8.9 ± 0.4 


11.4 ±0.4 


6.4 X 10- 


15 


0.27 


11917 ±402 


8.21 ± 0.08 


28 


H681 


232 ± 12 


4.36 ± 0.20 


260 ± 12 


9.7 ± 0.3 


27.7 ± 0.7 


2.6 X 10- 


14 


0.59 


14237 ± 559 


7.80 ± 0.09 



Note. — Line fluxes are normalized to H^=100, after correcting for reddening. F(H^) is the measured flux, collected for extinction. 
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Fig. 2. — Comparison of line fluxes and derived quantities for Hll regions from our new observations and those by |Kennicutt et ar]j2003| KBG03). (a) 
[OlI]A3727 (crosses) and [OlII]A5007 (asterisks); {h) [NlI]A6583 (diamonds) and [S II] AA67 17,6731 (triangles); fcj log(R23), {d) rj-feed oxygen abun- 
dances. The dashed line represents the one-to-one correlation. 
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TABLE 3 

Comparison with previous observations 









H972'' 


H237'' 


H149'' 


H140'' 


H681'' 




KBG03 


188 


148 


136 


212 


317 


244 


[Oll]A3727 


B07 


221 


139 












New 






192.8 


238.5 


351.8 


232 


[Oiii]A5007 


KBG03 


103 


30 


114 


318 


114 


304 




B07 


102 


30 












New 






106.1 


328.6 


145.4 


259 


[Nii]A6583 


KBG03 


64.6 


86.9 


48.3 


35.9 


60.2 


11.3 




B07 


71 


80 












New 






47.5 


35.6 


57 


9.7 


[Sii]AA6717,6731 


KBG03 


28.8 


38.4 


35.3 


37.2 


80.9 


27 




B07 


29.8 


41 












New 






38.5 


39.4 


67.1 


27.7 


[Oiii]A4363 


KBG03 








1.8 








B07 


0.24 














New 








1.7 






R23 


KBG03 


3.27 


1.88 


2.90 


6.40 


4.71 


6.54 




B07 


3.58 


1.79 












New 






3.36 


6.81 


5.48 


5.82 


12+log(0/H)j; 


KBG03 


8.71 






8.33 




7.92 


B07 


8.52 














New 








8.40 




7.80 



Note. — Six H II regions are in common between KBG03, B07 and our new observations. 

^ For H1013 and H972, the B07 line fltixes were adopted to construct the enlarged H 11 region sample 

(see text). 

For H237, H149, H140 and H681 the newly-determined Une fluxes from our observations were used 
for the enlarged sample. 
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Fig. 3. — Nebular diagnostic diagrams showing tlie excitation se- 
quence for our Hll region sample: log( [N II]A6583/Hq) (left) and 
log([SlI]AA6717,6731/Ha) (right) as a function of log( [O III]A5007/H/3). 
The curves re present upper bounda ries for star-forming r egions ionized by 
hot sta rs fromlKewley et al.H2001| long-dashed line) and |Kauffmann et al.| 
j2003| short-dashed line). The tilled blue triangles mark the positions of HI/ 
in our new observations. Symbols as in Fig.[T| 

H1013, B07 provided an oxygen abundance 0.21 dex lower 
than KBG03. We adopted the value from B07, obtained from 
deeper data, and which corresponds to an zone temper- 
ature directly derived from [O lll]A4363/(A4959 -I- A5007), 
while in KBG03 it was derived from T[Slll]. Finally, we 
adopted the H972 measurements from B07, while for the re- 
maining four H II regions included in Table [3] our new obser- 
vations were adopted. 

Our final sample comprises 79 H II regions, with 28 objects 
from our new observations, 20 from KBG03, 4 from B07, 
16 from KG96 and 11 from vZ98. In this merged sample, 
there are 29 measurements of the [Olll]A4363 auroral line 
(10 from our new observations, 17 from KBG03 and 2 from 
B07). In order to ensure consistency in the analysis, we re- 
calculated the Te and the direct oxygen abundances when the 
[Olll]A4363 line was measured. 

In Fig. [3] we examine the excitation properties 
of our Hll region sample by plotting the classic 
[Oiii]A5007/H/3 vs. [Nii]A6583/Ha and [ Oiii]A5007/H^ 
vs. [S II]A6717,6731/Hq! diagnostic diagrams (Baldwin et al. 

:BPT). 



M9811: 



We include in the plot the boundaries between different 
photoionization sources (star-forming regions vs. AGN) by 

twley et al. (2001 shown by the long-das hed Une) and 
auffmann et al.„2003) short-dashed line). [Kewley et al. 
(|2001|l used stellar population synthesis and photoionization 
If 



models to create the 'maximum starburst line' on the BPT di- 
agrams. A modification to this classification was provided by 
|Kauffmann et al . ( 2003 ) to include an empirical boundary line 
between pure star-forming galaxies and Seyfert-H II compos- 
ite objects. As shown in Fig.p] according to these criteria all 
the objects in our sample are located in the pure star-forming 
region of the diagram. 



3. RADIAL ABUNDANCE GRADIENT AND LOCAL AZIMUTHAL 
VARIATIONS 

The first measurements of a radial abundance gradient from 
H II regions in MlOl can be traced back to more than 40 years 
ago (Searle 1971; Smith 1975). Based on our new spectro- 
scopic observations, in combination with the work published 
by other authors, as presented in Sect. |2] we are in the po- 
sition to better assess the spatial distribution of the oxygen 
abundance in the disk of this galaxy. Given the relatively 
dense spectroscopic coverage in the western part of the disk 
we can investigate how the oxygen abundance varies along the 
arc-like region identified in Fig.fT] designated as arc A here- 
after, covering the restricted radial range R/Ro = 0.54- to 0.59, 
and extending in position angle between 230 and 290 degrees 
(^17 kpc in projected length), where we have a good amount 
of nebular spectra with [O III] A4363 detections, in order to 
test whether we could confirm the local azimuthal variation 
suggested by KG96. 

3.1. Abundance variations based on the direct method 

There are 29 Hll regions with [Olll]A4363 line measure- 
ments in our sample, covering the radial range from R/Rq = 
0.10 (H493) to R/R() = 1.25 (SDH323). The resulting radial 
oxygen abundance gradient is presented in Fig.|4] As the fig- 
ure shows, the data can be well fitted by a weighted linear 
least-squares fit (the Hnear Pearson correlation coefficient is 
-0.90): 



l2 + log(0/H)T^ = 8.73(±0.03) - 0.87(±0.04) R/Rq, (3) 

which, unsurprisingly, is consistent with previous re- 
sults by KBG03, 12 -I- log{0/H)T = 8.76(±0.06) - 
0.90(±0.08) R/Ro, and B07, 12 -I- log(C>///)r = 
8.75(±0.05)-0.90(±0.07) R/Rq. 

Let us examine the abundance distribution along arc A. As 
shown by the gray shadowed region in Fig. |4] this part of the 
galaxy contains objects whose O/H abundance ratios appear 
to cover a larger spread (at constant radius) than elsewhere. 
To clarify the situation, in Fig.|5]we zoom into this restricted 
radial range (top panel), and plot the oxygen abundance as 
a function of position angle (bottom panel). The red dotted 
lines mark the range in O/H corresponding to the change in 
galactocentric distance from 0.54 TJq to 0.59 Rq, as inferred 
from the gradient given by Eq[3] 

The Hll region H27 stands out with an abundance 
12H-log(0/H) = 7.96 ±0.08, while the corresponding radial fit 
value at its galactocentric distance is 8.23 ±0.05. The dif- 
ference between the measured and the fit value is 0.27 ± 
0.09 dex. Therefore, H27 is a 3a outlier from the linear fit. 
An additional object observed by KBG03, H128, is also an 
outlier relative to the radial fit, with an observed value of 
12H-log(0/H) = 8.45 ±0.04 dex and a fit value of 8.25 ±0.05 
dex (^3.1cr). All other objects are consistent with the radial fit 
value for their galactocentric distances, within 3x their mea- 
surement uncertainties. 

Like all other objects, H27 and H128 are located in the 
pure star-forming region of the BPT diagram (see Fig. |3]l 
and we find no evidence for the presence of other ionization 
mechanisms (e.g. shocks) that could explain its peculiar oxy- 
gen abundance for its galactocentric distance. For example, 
we find no significant enhancement in the strength of shock- 
sensitive lines^su£has[S II], [O l]. H27 lies within the region 
in which [Franchetti et"aL] ( |2012j ) made a spatially-detailed in- 
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Fig. 4. — Radial oxygen abundance gradient determined from direct electron temperature measurements. The filled circles represent our new observations. 
The H II regions from KBG03 and B07 are labeled by the diamond and the asterisk symbols, respectively. The solid line is the weighted linear least-squares fit 
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culated from Eq 



oxy, 

■i 



vestigation of supernova remnant (SNR) candidates in MlOl 
from archival Hubble Space Telescope images, and has been 
identified as a superbubble, rather than a SNR. 

Analyzing the scatter relative to the radial abundance gra- 
dient is one app roach to probe azimuthal variations of metals 
( BresoUn|201 1 1. The rms scatter around the least-squares fit 
for our full sarnple of H II regions (black solid line in Fig. Hll 
is 0. 1 1 dex. For the 1 1 objects contained within the grey shad- 
owed portion such scatter rises to 0.15 dex, while for the 18 
remaining objects it is 0.08 dex. 

A scatter of 0.15 dex potentially represents a detection of 
local inhomogeneities in the abundance distribution. In order 
to test the significance of this result, we calculated the errors 
on the scatter measured for these two data sets (inside and 



, following ajackknife 



outside the shado wed reg ion in Fig. Hll, 
procedure (Lupton 1993' p. 46). Simply put, jackknife resam- 
ples a sample of size A^, constructing subsamples of size 
A^- 1 by omitting an element of the sample in succession, and 
then estimates the variance of the scatter from 



N 

■E 

!=1 



-tN-l) 



(4) 



where tj^-i is the average of t over the N subsamples. 

By applying this technique, the resulting scatters are 
0.15 ±0.03 dex (objects within arc A) and 0.08 ±0.01 dex 
(objects outside arc A), respectively. Thus, the significance 
of the difference is only at the ^2a level. 

In summary, using the direct method on our enlarged sam- 
ple of 29 H II regions with [O III] A4363 detections, we derive 
a radial oxygen abundance gradient that is consistent with pre- 
vious results. We detect two outliers from the mean radial 
gradient, H27 and HI 28, both of which with a significance at 
the 3(7 level. No evidence of shocked gas are found for H27. 
The objects located along arc A display an abundance scatter 
of 0. 15 ± 0.03 dex, which represents a marginal detection of a 
local azimuthal variation. 

3.2. Abundance variations based on strong-line methods 

3.2.1. Comparison of different strong-line methods 

In addition to the [O lll]A4363-based method, in order to 
take advantage of the 79 objects in our enlarged sample of 
H II regions we also considered a number of strong-line meth- 
ods, in which the strength of easily observed nebular emission 
lines can be used to calculate the nebular abundances. Among 
the available methods making use primarily of O lines we 
considered the R23 caUbration by McGaugh j (|1991| =M91 ) 
(in the analytical form provided by Kobulnic ky et al.|1999 l, 
and th e P-method by |Pilyugin| ( [^01, = POl ) and 'Pilyuginj 
( 2005 1 =P05). Other methods rely on the calibration of the 
strength of other emiss ion lines, typi cally the N lines, in terms 
of the O abundance: IPettini & Pagel (2004 PP04-N2 and 
PP04-O3N2), Perez-Montero & Contini ( 2009 PM09-N2 and 



PM09 -O3N2), |Pilyugin et al.,(,2010, ONS) „Kewley & Dopita 



( |20()2l KD02) and th e empIricaTcahbration of [N iij/ [O iij 
irom'Breso lin| ( |20(y7l NO). These methods are summarized in 



Table |4] for clarity. They are somewhat arbitrarily divided in 
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'empirically' and 'theoretically' calibrated methods, to reflect 
the fact that the calibration in terms of 0/H abundance ratio 
can be obtained from samples of H II regions with available 
direct abundances or from grids of photoionization models, 
respectively. 

For our purposes, we were interested to test how the strong- 
line abundances compared with the abundances from the di- 
rect method. We show in Fig. l6] the comparison between H II 
region abundances obtained from the Tg method (along the 
X-axis) and those obtained from the various strong-line diag- 
nostics. Since similar comparisons have appeared frequently 
in the literature in the past decade, presenting comparable re- 
sults, we keep our discussion limited to the essenti al. For 
fu rther d etails, the reader is re ferred to, fo r example: |Kenni^ 

flFeTal. (2003 ), Bresolin et al. (2009a^, 'Kewley & Eillson 
008) , Perez -Montero & Diaz (2005 ), and Yin et al. (2007) . 

Empirically calibrated methods - The first two rows of Fig. [6] 
refer to empirically calibrated diagnostics. A quick look at 
the figure shows that despite the fact that, by construction, 
these methods should generally agree with the direct method, 
the level of agreement varies greatly, at least for the sample 
of H II regions we are considering. In some cases this is due 
to the fact that the O/H range in which the diagnostics were 
calibrated is narrower than the full range of abundances avail- 
able in MlOl (which is unusually large). Secondly, possible 
differences in physical conditions (e.g. ionization parameter) 
between the calibrating sample and our sample can translate 



into the variations we observe. With these considerations in 
mind, we point out the following points: 

• PP04-N2 provides a relatively good match to the di- 
rect abundances (with a rms scatter of 0.11 dex), 
but large discrepancies arise at low metallicity, 
12H-log(0/H) < 8.0. Similarly for PP04-O3N2, but the 
scatter is worse (0.17 dex). 

• PM09-N2 tends to systematically underestimate the 
abundances, compared to the direct method. PM09- 
03N2 yields similar results. 

• the POl calibration was derived from nebular data with 
12-i-log(0/H) > 8.2, which is therefore the applicable 
range of this method. Accordingly, it is not surprising 
to see large discrepancies in the low metallicity regime. 

• P05 is a revised version of POl, providing calibrations 
for both the upper and the lower branch of the relation 
between oxygen abundance and Rn- Here, we only 
adopted the calibration for the upper branch, but even 
at high metallicity the diagnostic underestimates the T^- 
based abundances. 

• the ONS calibration provides abundances that are com- 
parable to those from the method across the whole 
metallicity range, with a rms scatter of 0.09 dex. 
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• the NO calibration seems to largely agree with the 
method, although the scatter of the data points is larger 
than in the case of the ONS method. 



Theoretically calibrated methods - The last row of Fig. [6] 
presents the abundances obtained from strong-line methods 
that are calibrated from theoretical models. It is well-known 
that they provi de systematically high er abundances than the 
Te method (e.g. |Bresolin et al.|2009 aV 'KD02-combine' rep- 
resents the so-called 'combined' diagnostic of KD02, recom- 
mended as an optimized abundance determination by these 
authors, and which can be used over the range of abundances 
from 12-i-log(0/H) = 8.2 up to 9.4, with a minimized scat- 
ter and showing no systematic offset compared with 'KD02- 
average'. 'KDO2-R23' is the R23 method calibrated from the 
models of KD02: compared to the method it shows a large 
scatter both at high and low metallicities. 'KD02-average' is 
the average of five independent strong-line methods. As we 
can see in Fig. [6] it overestimates the T^-based abundances, 
but the relative onset appears unchanged across the full metal- 
licity range. 

Finally, M91 suffers from the well-known double-valued 
nature of the R23 diagnostic. Additional line ratios and/or 
an initial guess of the abundance are required to choose the 
appropriate branch. We attempted to break the R23 degen- 
eracy by using the [N ll]/ [ O ll] and [N Il]/Ha ratio s, fol- 
lowing the criteria found in |Kewley & Ellis"on ( 2008| l. Ac- 
cording to these authors, the separation between the R23 up- 
per and lower branches occurs at log( [N ll]/ [O ll]) = -1 .2 and 
-1.3 <log( [Nll]/Ha)<-l.l. These boundary values are 
marked by dashed lines in Fig. 17] where we plot the trends 
of both [N II]/ [O II] and [N Il]/Ha with R23. Only when both 
these criteria were simultaneously satisfied we assigned the 
objects to the appropriate branch. 



However, there is a large number of 22 objects in our sam- 
ple whose branches could not be firmly established, or remain 
dubious or contradictory, based on these criteria. These ob- 
jects are represented by black filled circles connected by dot- 
ted lines in Fig. [8] where we show how the M91 -derived abun- 
dances change as a function of R23 parameter (top panel) and 
galactocentric distance (bottom panel). As the bottom panel 
of this figure shows dramatically, the R23 branch choice can 
yield very different conclusions about the radial trend of the 
oxygen abundance. This can make double-valued diagnos- 
tics like R23 inappropriate to study radial abundance gradients 
in galaxies, and as a minimum the results should be checked 
with additional diagnostics. 

In summary, the ONS calibration gives the best agree- 
ment with the Te method for the oxygen abundances of the 
MlOl Hll region sample among the strong-line diagnostics 
considered in our study. Most of the remaining empirically- 
calibrated strong-line methods diverge from the method 
at low metallicities. The theoretically-calibrated diagnostics 
tend to systematically overestimate the abundances, com- 
pared to the method. The excellent agreement of the ONS 
method with the method is not too surprising, given that 
it was specifically calibrated using Hll regions with well- 
measured Tg-based abundances (but so where the others). To 
further test its agreement with the method we applied the 
ONS method to the H II reg ion sample in the galaxy NGC 300 
presented by Bresolin et al. (2009a), who found good agree- 
ment between the nebular abundances derived with the 71 



method and the m etallicities of blue supergiant stars ( Ku-| 
|dritzki et al.||2008y The abundance gradient from the ONS 
method is found to be represented by the following linear 
fit: 12H-log(O/H) = 8.49(±0.02) - 0.33(±0.04) R/Rq. 
This agrees well with the fit obtained from the Tg 
method, 12H-log(0/H) = 8.57(±0.02) - 0.41(±0.03) R/Rq, 
and from the blue supergiants, 

12H-log(O/H) = 8.59(±0.05) - 0.43(±0.06) R/Rq. 

3.2.2. Abundance variations based on the ONS diagnostic 

In order to expand our analysis of possible spa tial v ariations 
of the chemical abundances presented in Sect. |3.1| we used 
the oxygen abundances derived from the ONS niethod, which 
shows the best agreement with the Tg method, for the enlarged 
sample of Hll regions. However, adopting other strong-line 
diagnostics would yield the same conclusions concerning az- 
imuthal variations in the inner disk. 

The radial oxygen abundance gradient thus derived is 
shown in Fig. [9] The solid line represents the least-square fit 
to the data. We note that there are five objects at /?//?() = 0.81 
displaying a large abundance scatter. These are five individ- 
ual emission knots located within the supergiant H II region 
NGC 5471 (see KG96). At least two of them have been 
associ ated with SNRs or the presence of high-velocity gas 
(Skill man||1985j |Chu & Kennicuttl[T986l [Kennicutt &"GEn 
nett 1996|l. Peculiar ionization conditions in these knots could 
in principle affect the chemical abundance determination (but 
we note that the Tg-based oxygen abundances do not appear to 
be affected, and are quite homogeneous between the different 
knots). Thus, the apparently large abundance scatter between 
these data points, as seen in Fig. [9] is unlikely to result from a 
real oxygen abundance inhomogeneity in this region (an inter- 
pretation supported by the Tf-hased O/H measurements; see 
KBG03 for a discussion of the N/O ratio), but rather from a 
breakdown of the ONS diagnostic in this nebular complex. 
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TABLE 4 

Summary of nebular abundance diagnostics adopted 



Number 


Method 


Emission lines involved 


Calibration class 


1 


Te 


[Olll]A4363, [Olll]A4959,5007 


Direct 


L 




[IN llJ/rlQ; 


Empirical 


3 


PM09-N2 


[Nll]/Ha, [Nil]/ [On] 


Empirical 


4 


PP04-O3N2 


[Nll]/Ha, [OlII]/H/3 


Empirical 


5 


PM09-O3N2 


[Nll]/Ha, [OlII]/H/3,[NlI]/[OlI] 


Empirical 


6 


POl" 


[O III]/ [Oil] 


Empirical 


7 


P05" 


[O III]/ [On] 


Empirical 


8 


ONS 


[S II]/ [Oil], [Nil]/ [Oil], [OlII]/H/3 


Empirical 


9 


NO 


[NlI]/[0 II] 


Empirical 


10 


M91 


R23, [O III]/ [Oil] 


Theoretical 


11 


KD02-combine'' 


[NlI]/[OlI], R23, [OlII]/[OlI], [NII]/[SII] 


Theoretical 


12 


KDO2-R23" 


[NlI]/[OlI],R23, [O III]/ [On] 


Theoretical 


13 


KD02-average'' 


[NlI]/[OlI],R23, [O III]/ [On] 


Theoretical 
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Fig. 8. — Top: oxygen abundances from R23 as a function of log(i?23). Bottom: radial oxygen abundance distribution based on the M91 calibration of the S23 
diagnostic. There are 22 objects (black filled circles) whose branches could not be univocally established and both upper and lower branch O/H values are shown 
as black dots connected with dot lines. Symbols as in Fig.[T] 



The linear fit to the radial abundance gradient is not signif- 
icantly affected by the inclusion of the knots in NGC 547 1 . 
We obtain: 



12 + log(0///)oiV5 = 8.74(±0.03)-0.82(±0.05) R/Ro (5) 

■with a linear Pearson correlation coefficient of -0.90. The 
intercept and the slope of the regression are consistent ■with 
the result derived from the direct method (Eq.[3|. 

As done earlier, in Fig.[TO](top panel) ■we zoom into the ra- 
dial range R/Rq = 0.54 to 0759 (highlighted in gray in Fig.|9]l, 
and plot the oxygen abundance as a function of position an- 
gle (bottom panel). Hll regions H27 (filled blue triangle) 
and H128 (filled green circle), ■which were found to be possi- 
ble outliers from the radial abundance gradient obtained from 



the direct method, are not so using the ONS method. Ho^w- 
ever, a new object, corresponding to slit 6 in |van Zee et al. 
( |1998 1, is now found to have quite a low O/H abundance ra- 
tio, 12-i-log(0/H) = 7.84. This can be explained by the fact 
that the criterion used in the ONS method to assign an H 11 
region to a particular excitation class fails for this particu- 
lar nebula. The criterion, based on the [N Il]/[S ll] line ra- 
tio, assigns this object to the 'hot' class, but if we calculate 
its O/H abundance as if it were a 'warm' H II region, its O/H 
ratio would increase by 0.3 dex, which would thus remove 
the systematic abundance offset relative to the other nebulae 
with similar galactocentric radius. The problem of misclas- 
sification in the ONS method has already been discussed by 
IB resoUn et al. ( 2012 1 . In a ddition, for the specific case of 
MlOl Pilyugin et al. (2010 1 recommended to use the "hot' 
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Fig. 9. — Radial oxygen abundance gradient from tlie ONS method. The solid 
the filled triangle. The grey shadowed region is the same as in Fig. [4] Symbols as 



o 

o 



line is the linear least-squares fit including 74 H II regions. H27 is marked by 
in Fig.[T] 
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Fig. 10. — Detailed view of the shadowed region in Fig.|9] The upper panel 
shows the relation between oxygen abundance and galactocentric distance in 
the range R/Ro = 0.54 to 0.59. The bottom panel is the relation of oxygen 
abundance with position angle. The filled circle and triangle represent H128 
and H27, respectively. Symbols as in Fig.[T] 

classification only for objects lying at R/Ro > 0.90. 

The rms scatter around the least-square fit for 74 H II re- 
gions (the enlarged sample with the five knots in NGC 5471 
removed) is 0. 10 dex. For the 17 objects in the grey shadowed 
region of Fig. j9]the scatter is 0. 14 dex, and for all the remain- 
ing objects it IS 0.09 dex. However, if we focus on the 13 ob- 
jects along arc A only we find that the scatter is comparable to 
what we find for the remaining objects (^0.10 dex). In other 
words, based on the ONS method, we find no evidence for 
local inhomogeneities of the oxygen abundance along arc A. 

4. A GLOBAL ASYMMETRY IN THE OXYGEN ABUNDANCE? 

It is interesting to test not only for local chemical inhomo- 
geneities, as done in the previous sections, but also for large- 
scale azimuthal variations in the disk of MlOl. Is it true that 
the metal abundance distribution in the disk is azimuthally ax- 
isymmetric, as usually assumed in chemical evolution mod- 
els? To start to answer this question in the case of MlOl, 
further examination of the two dimensional distribution of the 
oxygen abundance is required. 

As a simple test, we divided the galaxy into two halves with 
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Fig. 1 1 . — Spatial distribution of the H II regions in our sample in the plane 
of the sky. Our sample is geometrically divided into the 'SE' (with position 
angle relative to the galaxy center between 37 and 217 deg) and 'NW part 
(complementary values of position angle). Symbols as in Fig.[T] 
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Fig. 12. — Radial dependence of the abundance parameter R23- Filled 
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respect to the major axis: a SE part (position angle: 37 to 217 
degrees, including 37 objects) and a NW part (complemen 
tary position angles, with 42 objects), as shown in Fig. fTT 



KG96 suggested that objects in the SE have lower oxygen 
abundances compared to H II region in the NW. This result 
was obtained from the use of the R23 parameter. We plotted 
the radial dependence of R23 in MlOl in Fig. 12 in the same 



way KG96 did, but including about twice the number of ob- 
jects (79 versus 37). The distribution of data points is virtually 
the same for the two subsamples, and thus we do not confirm 
the non-axisymmetric distribution suggested by KG96 con- 
sidering the /?23 diagnostic alone. 

As a further test, we used six different strong-line 
methods to calculate the oxygen abundances and fitted 
the radial distribution with linear functions of the form 
+ bx{R/Ro) 
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12H-log(0/H) = fl + bx(R/R^i). Since we find (see Fig 
that different abundance diagnostics yield qualitatively differ 
ent radial trends in the outer disk, only objects with 7?/7?()<0.8 
were included in the fit. In Table |5] we summarize, for each 
diagnostic, the linear fit parameters a and b and their uncer- 
tainties, the Pearson correlation coefficients, the rms scatter of 
the fit and the characteristic oxygen abundances measured at 
R = OAR(), which correlate with the integrated galaxy metal- 

licity (Zaritsky et al. 1994, Moustakas & Kennicutt 2006jl. galacto centric distance, as already found by Kennicutt et al. 



the presence of azimuthal variations among smaller sections 
(e.g. quadrants) of the disk. 

5. SUMMARY AND CONCLUSIONS 

Using a data sample of 79 Hll regions, with 28 from 
our new observations (yielding 10 new detections of the 
[Olll]A4363 line) and the rest from additional sources in 
the literature, we have obtained and analyzed the radial 
oxygen abundance gradient in MlOl based on the direct 
method ( [O ni]A4363-based), as well as various strong- 
line diagnostics. We found an exponential abundance pro- 
file with a gradient of -0.87 ±0.04 dex Rq"' and a cen- 
tral abundance of 12H-log(0/H)=8.73 ±0.03, from the direct 
oxygen abundance measurements. The scatter in the ra- 
dial abundance gradient along arc A, located in the west- 
ern disk of MlOl, is 0.15 ±0.03 dex, while for the remain- 
ing Hll regions of the MlOl disk the scatter is 0.08 ±0.01 
dex. In the same section of the galaxy we found that one 
Hll region, H27, deviates from the overall galactic abun- 
dance gradient, having a significantly lower oxygen abun- 
dance, 12H-log(O/H)j-e=7.96±0.08, compared to nearby ob- 
jects. One additional neighboring nebula, H128, has instead 
a large metalHcity value, 12H-log(O/H)7-e=8.45±0.04, for its 



Fig. 13 illustrates the results of the fitting procedure for the 



two hTi region subsamples (SE and NW) for each diagnos- 
tic. Table [5] shows that the linear radial gradient parame- 
ters (e.g. slope and zero point) are highly dependent on the 
abundance diagnostic adopted, which is not surprising given 
the different behaviors found among the strong-line methods 
when comparing with the method. The possible breakdown 
of many of the diagnostics at low metallicities is evident in 
Fig. 



( 2003 1. These results provide evidence that marginally signif 
icant deviations from local metallicity homogeneity can arise 
in the interstellar medium of this part of the galaxy over spa- 
tial scales of a few kpc. 

Among the strong-line abundance diagnostics we consid- 
ered to derive the O/H abundance ratio for the full sample 
of 79 H II regi ons we found that the ONS method (Pilyi>] 
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which shows that in several cases (e.g. for POl and 
KD02^ombine) the radial abundance gradient reaches a min- 
imum abundance at ^0.8 /?o, and then rises again in the outer 
disk. This situation has been more thoroughly discussed by 



gin et al. 2010 1 provides the best agreement with the direct 



Pilyugin ( p003j ). If we believe that these results are unphys- 
ical, i.e. that the actual galactic gradient does not 'turn over' 
at a certain galactocentric distance, we need to discard the 
abundances inferred from the strong-line diagnostics that are 
affected by this issue (we stress that the direct abundance are 
not affected). However, we point out that in recent years sev- 
eral reports of metallicity gradients becoming flat or turning 
over at large galactocentric radii in the Milky Way and other 



galax i es have been made (e.g . Worthey et al. 2005 ; Yong et al. 
|2012[[Bresolin et al.|2009bl l. The radius where such a break 
appears to occur in MlOl agrees well with the results from the 
recent work of Scarano & Lepine (2013), who noted a good 
correlation between the radial positions of the break and coro- 
tation radii for a sample of spiral galaxies. Our Fig. 13 shows, 
however, that the various abundance diagnostics we consid- 
ered still leave some ambiguity concerning both the presence 
of the break and the gradient slope in the outer disk (this prob- 
lem is absent in other outer disk nebular abundance studies. 



e.g. |Bresolin et a l. 2009b, 20121. Since the main goal of this 
work is to investigate the possibility of azimuthal variations in 
the inner disk, where the metallicity trends are less ambigu- 
ous, we decided to limit our analysis to R/Rq < 0.8. 

Doing so, the differences in slope and intercept of the lin- 
ear regression between the SE and the NW portions of the 
disk are found to be smaller than the corresponding uncertain- 
ties, confirming the absence of obvious large-scale abundance 
variations in the azimuthal direction. A much more extensive 
spectroscopic coverage of MlOl would be required to test for 



method. With the ONS method we did not find any signifi- 
cant difference in the abundance scatter between H II regions 
in the western arc of MlOl and the rest of the galaxy, being 
approximately 0.10 dex in both cases. The rather large devi- 
ations from the main abundance radial gradient measured for 
H27 and HI 28 from the direct method are not confirmed by 
the ONS method. On the other hand, the ONS method pro- 
vides a large (0.4 dex) spread in the abundance of individual 
knots in the supergiant Hll region NGC 5471, which is not 
seen in the T^-hased data, and which we interpret as an exam- 
ple of the fact that strong-line method results should always 
be taken with care, and be considered valid in a 'statistical' 
sense. For individual measurements the Tg method should be 
preferred. On the other hand, the result of our test of possible 
large-scale azimuthal abundance trends between two opposite 
sides of MlOl using strong-line methods, i.e. that there is no 
detectable variation, should be considered robust. The small 
scatter we generally observe in the abundance gradient, when 
we consider that it includes H II regions on opposite sides of 
the galaxy, is also a good indication for a fairly homogeneous 
azimuthal abundance distribution. 

In general, findings of large local deviations from abun- 
dance homogeneity obtained from strong-line methods should 
be verified with direct, 7^.-based measurements. We point out 
that the abundances of individual emission knots within large 
(^1 kpc) Hll region complexes are consistent with chemi- 
cal abundance homogeneity when using r^-based data. Be- 
sides the case of NGC 547 1 mentioned above, this is true for 
the NGC 5447 (H128, H143, H149) and NGC 5462 (HI 159, 
HI 170, HI 176) complexes. On the other hand, claiming siz- 
able inhomogeneities from abundances measured from some 
strong-line indicators appears risky. For example, the N2 pa- 
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TABLE 5 

Parameters of the separate linear fits to the radial gradients 



Region 



(T(a) a{b) correlation coeff. scatter (rms) 12+log(0/H)K=o 



PP04-N2 

ALL 

SE 

NW 



!.72 
!.66 



PP04-O3N2 

8.85 



ALL 
SE 
NW 

POl 

ALL 

SE 

NW 

ONS 

ALL 

SE 

NW 

NO 

ALL 

SE 

NW 



!.84 



8.77 
8.72 
8.82 



8.75 
8.74 
8.77 



8.75 
8.73 
8.78 



KD02-combine 

ALL 9.15 
SE 9.11 
NW 9.19 



-0.75 
-0.88 
-0.69 



-1.03 
-1.15 
-1.00 



-0.78 
-0.69 
-0.86 



-0.85 
-0.88 
-0.87 



-1.18 
-1.22 
-1.22 



-0.96 
-0.90 
-1.02 



0.03 
0.04 
0.04 



0.05 
0.06 
0.07 



0.03 
0.04 
0.05 



0.03 
0.05 
0.04 



0.03 
0.04 
0.04 



0.03 
0.04 
0.04 



0.06 
0.10 
0.09 



0.09 
0.15 
0.13 



0.07 
0.10 
0.09 



0.06 
0.11 
0.07 



0.06 
0.10 
0.08 



0.06 
0.10 
0.07 



-0.83 
-0.87 
-0.79 



-0.80 
-0.84 
-0.77 



-0.81 
-0.79 
-0.82 



-0.87 
-0.84 
-0.89 



-0.92 
-0.93 
-0.93 



-0.90 
-0.88 
-0.92 



0.10 
0.08 
0.10 



0.14 
0.13 
0.15 



0.10 
0.09 
0.11 



0.09 
0.10 
0.09 



0.09 
0.08 
0.09 



0.08 
0.08 
0.08 



8.38 
8.37 
8.38 



8.44 
8.42 
8.45 



8.46 
8.45 
8.48 



8.41 
8.39 
8.42 



8.27 
8.25 
8.30 



8.77 
8.75 
8.78 



Note. — Linear fits a + fcx(R/So) including only objects withR/So < 0.8 



rameter (in both the Pettini & Pagel[2004 and Perez-Montero 
|& Contini||2009 1 versions of the calibration considered here) 
yields abundance variations on the order of ~0.2 dex between 
the few knots contained in the Hll region complexes men- 
tioned above, i.e. over spatial scales of a few hundred pc. 
Among the N-based diagnostics, the NO indicator calibrated 
by Bresolin ( 2007 ) appears much more robust, yielding vari- 
ations on the order of 0.03 dex among these objects. 

It is tempting to try and interpret the relatively small abun- 
dance peculiarity we detected in the western disk of MlOl 
using the direct method with the notion that metallicity inho- 
mogeneities could arise from metal-poor gas infall and galaxy 
interactions. The presence of high-velocity clouds and distor- 



sult about the errors on the scatters, R.P. Kudritzki for useful 
comments, and T.-T. Yuan for her help with various scien- 
tific discussions. FB gratefully acknowledges partial support 
from the National Science Foundation grants AST-0707911 
and AST- 1008798. We thank the anonymous referee for con- 
structive comments that helped us to improve the quality of 
the paper. 



tions in the spiral structure in the eastern disk of MlOl (van 
der Hulst & Sancisi|[T988l [Waller et al.|[T997l ISancisiltaT 



2008| l can be attributed to tidal interactions between MlOl 



and one or more o f its surround i ng nei ghboring dwarf galax- 
ies. Very recently, |Mihos et al. ( |2012p presented a new deep 
Hi map of the MlOl group, which revealed the presence of 
newly-discovered Hi clouds and an Hi 'plume' extending 
from the SW of MlOl, likely of tidal nature. It is of course 
difficult to associate any of these features with the observed 
metallicity patterns in the MlOl disk, but it is tantalizing to 
observe that the marginally larger scatter in the abundance 
distribution that we detect in the western arc is taking place 
on the same side of the galaxy where these tidal features ap- 
pear now more evident. Future studies should address this 
possibility by obtaining a large number of high-quality spec- 
troscopic observations and direct abundances in this part of 
MlOl. 
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Fig. 13. — Radial oxygen abundance gradients in the SE and NW sections of MlOl based on six different strong-line methods. The straight lines show the 
weighted linear regressions to the two samples, limited to data points located at R < 0.8 Rq (marked by the vertical dashed line). 
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